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Abstract. 

Observations of the nebula S3 associated with the WO star in the galaxy IC 1613 and of an 
extended region surrounding S3 are reported. The star and bright core of the nebula were observed 
with a multipupil fiber spectrograph mounted on the 6-m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences. Images in the principle spectral lines and 
integrated spectra of the star and three compact clumps were obtained, and the radial-velocity 
field constructed. An extended region of the galaxy was observed with the Very Large Array at 
21 cm. A giant ring or H I shell enclosing a large fraction of the stellar population in IC 1613 was 
discovered. The WO star and associated bipolar nebula, which we discovered earlier, lies at the 
inner edge of the H I ring. A local H I deficiency and two arc-like H I ridges were also detected 
for the first time, and probably represent the neutral component of the bipolar shell surrounding 
the WO star. The two arc-like ridges may also have been produced by the collective stellar wind 
(and supernova explosions?) in OB association No. 9 from the list of Hodge. A scenario for the 
formation of the extended bipolar feature is discussed, based on the new data. 



1. Introduction 

Only six WR stars of more than five hundred such 
objects in the Local Group of galaxies belong to the 
rare WO class, which represents a very brief final 
evolutionary stage of massive stars approaching the 
naked CO-core stage (see Barlow & Hummer (1982), 
Kingsburg et al. (1995) and references therein). WO 
stars are characterized by powerful stellar winds with 
terminal velocities reaching ~ 5000 kms~^ ( Barlow 
& Hummer, 1982; Torres et al., 1986; Dopita et al., 
1990; Potcaro et al., 1992 ) and their effective tem- 
peratures can be as high as 10^ K (Dopita et al, 
1990; Maeder & Meynet, 1989; Melnik & Heydari- 
Malayeri, 1991; Polcaro et al., 1991). One of these 
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six stars was discovered by D'Odorico & Rosa (1982) 
(see also Davidson & Kinman 1982; Armandroff & 
Massey, 1985) in the dwarf irregular galaxy IC 1613. 
The star was identified by the specific wide bands it 
produces in the spectrum of the core of the bright 
emission nebula S3 from the list of Sandage (1971) - 
a giant 29" x 9" (93 x 30 pc) H II region. Throughout 
this paper, we adopt the distance 660 kpc for IC 1613, 
as derived by Saha et al. (1992). The central region of 
S3 shows bright He II 4686A emission (D'Odorico & 
Rosa, 1982; Smith (1995), Garnett et al., 1991). A de- 
tailed spectral analysis of the nebula and a discussion 
of its chemical composition are presented by Garnett 
et al.(1991) and Kingsburgh & Barlow (1995) - here- 
after KB95. The nebula S3 was found to be a source 
of thermal radio emission (Goss & Lozinskaya, 1995); 
the full size at half-maximum intensity of the radio 
source (19" x 14") corresponds to a bright region of 
the optical nebula visible on deep Ha images taken 
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Table 1: Log of MPFS observations. 



Date 


Size of 
micropupil 


Field of view 


Wavclcnth 
interval 


5X 


Exposure 
s 


Seeing 


29.XI.99 


1" 


16" X 15" 


4450 - 7050A 


8A 


7200 


2."7 


30.XI.99 


0."75 


12" X 11.3" 


4450 - 7050A 


81 


7200 


l."8 



by Hodge et al. (1990) and Hunter et al. (1993). 

Our observations made with a scanning Fabry- 
Perot interferometer mounted on the 6-m telescope 
of the Special Astrophysical Observatory (SAO) and 
narrow-band Ha and [O III] images taken with the 
4-m telescope of the Kitt Peak National Observa- 
tory resulted in the discovery of a faint external 
bipolar feature (hereafter called the bipolar shell) 
extending far beyond the bright nebula S3 (Lozin- 
skaya, 1997; Afanasiev et al., 2000 - Paper I). The 
sizes of the two external components^he southeast- 
ern and northwestern shells arc 35" x 24" (112 x 77 pc) 
and ~ 60" x 70"(190 x 220 pc), respectively Radial- 
velocity measurements at the maximum and the Ha 
line halfwidths provide evidence for expansion of the 
bipolar shell, with lower limits to the expansion ve- 
locities of the soiitheastern and northwestern shells 
being 50 and at least 7Gkms^^, respectively. 

In Paper I we proposed a scenario for the forma- 
tion of the unique bipolar feature by the stellar wind 
of the WO star, which is embedded in a dense layer at 
the boundary of a supercavity in the galaxy's H I dis- 
tribution. We identified this supercavity by analyzing 
the H I column density distribution in IC 1613 pub- 
lished by Lake and Skillman (1989). At that time, we 
hypothesized the presence of a dense shell surround- 
ing the supercavity and a physical association with 
the WO star, but were not able to prove that this 
was the case. 

In the current paper, we report the results of new 
observations of the nebula associated with the WO 
star in the galaxy IC 1613. We observed the bright 
core of the nebula S3 using the MPFS multipupil 
fiber spectrograph mounted on the SAO 6-m tele- 
scope (Section 2). We obtained integrated spectra of 
individual compact clumps in the central bright re- 
gion, and reconstructed the radial- velocity field from 
our H/3 and [O HI] 5007A measurements. We also 
obtained observations of a large sector of IC 1613 at 
21 cm with the Very Large Array (VLA; Section 3). 
We derived a detailed image of the neutral-gas dis- 
tribution, which revealed a giant dense ring or H I 
shell surrounding the supercavity, and showed that 
the WO star does indeed reside in this dense ring and 



is physically interacting with it. We also found a local 

H I deficiency and two extended features, testifying 
to the presence of a neutral gas component associ- 
ated with the ionized bipolar shell. In Section 4, we 
disciiss the results of the observations and analyze a 
possible scenario for the formation of the extended 
bipolar feature. 

2. Spectroscopic observations of the 
WO steir and nebula S3 

2.1. Observations and Data Reduction 

The spectroscopic observations were made on Novem- 
ber 29 and 30, 1999 using the new integral-field 
spectrograph MPFS (multipupil fiber spectrograph) 
mounted at the primary focus of the SAO 6-m tele- 
scope. A description of the spectrograph can be found 
at the web-page 

http : //www. sao . ru/~gaf an/devices .htm 

Compared to the earlier version of the Afanasiev 
et al.(1990), the new spectrograph has a larger field 
of view, wider spectral range, and higher quantum ef- 
ficiency. The spectrograph uses a 1024 x 1024 CCD 
as a detector and allows a set of spectra to be taken 
simultaneously from 240 spatial elements ('microp- 
upils') forming a 16 x 15 array in the plane of the sky. 
The angular size of an image element can be varied 
from C'S to 1.0". A spectrum of the sky background 
4f5 from the center of the field of view is taken simul- 
taneously. Table 1 gives a log of the MPFS observa- 
tions. 

These observations were reduced using IDL-based 
software developed by the SAO Laboratory for Spec- 
troscopy and Photometry of Extragalactic Objects. 
The preliminary data reduction included bias-frame 
subtraction, fiat-field correction, cosmic-hits removal, 
extraction of individual spectra from the CCD im- 
ages, and wavelength calibration using the spectrum 
of a He-Ne-Ar lamp. We then subtracted the night- 
sky spectrum from the linearized spectra and used a 
number of spectrophotometric standards observed at 
the same zenith angle as IC 1613 to convert the ob- 
served fluxes to an absolute energy scale. We used the 
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same standards to derive the seeing estimates given 
in Table 1. 

Analysis of the observations showed that the data 
obtained on November 29 and 30, 1999 yielded con- 
sistent values for the main parameters measured (line 
fluxes and velocities). Therefore, we will consider fur- 
ther only the November 30, 1999 data, since they had 
higher angular resolution. 

Figure la shows an [O III] line image of the central 
region of the nebula S3 surrounding the WO star from 
the data of Paper I, together with the boundaries of 
the field of view of our new MPFS observations. 

2.2. The WO star. 

We obtained the integrated spectrum of the WO star 
by co-adding 24 individual spectra, which exhibited 
the broad emission bands characteristic of Wolf-Rayet 
stars. 

During the flux measurements, we fitted the emis- 
sion lines using superpositions of two to six Gaus- 
sians; the resulting flux and equivalent- width esti- 
mates are summarized in Table 2. Note that the errors 
quoted in this table do not include uncertainty in the 
continuum level, which could contribute systematic 
errors to the equivalent widths. For comparison, this 
table also gives the fluxes and equivalent widths of 
the corresponding lines from KB95. 

We observed a narrower spectral interval than 
KB96; it is clear from Table 2 that the two sets of ob- 
servations do not show significant discrepancies. The 
differences between the equivalent widths are due pri- 
marily to the different methods used to determine the 
continuum level (unlike KB95, we did not take into 
account the nebular continuum component). Our fiux 
estimates are higher than those of KB95 for most 
of the lines. This is likely due to the fact that our 
panoramic detector collected all the light from the 
star, whereas a long slit whose width (1") was com- 
parable to or smaller than the image size was used 
in KB95. To test this hypothesis, we integrated the 
spectra of the star only over the lenses lying within 
a 1" X 5" rectangle aligned along PA = 55°, i.e., in 
the direction of the slit used in KB95. The resulting 
WO-star line fluxes, indeed, proved to be a factor of 
1.5-2 lower, on average, than the fluxes listed in Ta- 
ble 2. Moreover, as we show below, the emission-line 
brightness peak in the nebula does not coincide with 
the position of the star. This could have led to an 
offset during the positioning of the narrow slit on the 
star in KB95, resulting in further underestimation of 
the line fluxes. 

We constructed images in the C IV-|-He II 4659A, 



O V 5590 A, and C IV 5801 A emission hues by inte- 
grating the spectra over the intervals 4600 — 4725, 
5535 — 5653, and 5735 — 5850A, after subtraction 
of narrow nebular line components blended with the 
WO-star lines. The centroids of all the images agree 
to within ±0.05", and also agree with the ccntroid of 
a continuum image made in the vicinity of 6IOOA. We 
estimate the FWHM of these emission-line images of 
the WO star to be 1.85 ± 0.05", 1.84 ± 0.09", and 
1.67 ±0.03", respectively, consistent with the seeing 
estimate of 1^8 obtained for the standard star (Ta- 
ble 1). Note that, due to certain specific features of 
the spectrograph, the images constructed from the 
micropupil array are slightly elongated (1:1.15) along 
the Y axis (Fig. lb). Therefore, all our FWHM esti- 
mates are based on measurements in the X direction. 

The continuum image of the WO star was some- 
what more extended (FWHM = 2.2±0.15") and ex- 
hibited a faint base. (This could indicate the presence 
of a compact cluster made up of several stars, but firm 
conclusions require further observations). 

2.3. Structure and emission spec- 
trum of S3 

We constructed an entire series of monochromatic 
images of the bright core of the nebula in the 
He II 4686A, [O HI] 5007A, [S II] 6717±673lA, and 
Ha lines by obtaining Gaussian fits to the emission 
lines. The results arc shown in Figs. Ic-lf. All the 
line profiles are fit well by single Gaussians, and 
do not show any systematic deviations. Note, how- 
ever, that the expected relative gas velocities (up to 
100 km s^^ in the bright core Paper I) are apprecia- 
bly lower than the spectral resolution of our data 
(FWHM = 350 - 450kms-i). 

It is evident from Fig. 1 that the morphology of 
the nebula is appreciably different in different lines. 
For example, the He II brightness peak is coincident 
with the position of the WO star; the He II brightness 
decreases with distance from the star, and the full 
size of the emission region is 10" x 9". The situation 
is quite different in the brighter [O HI] and Ha lines, 
and also in the [S II] lines: the corresponding images 
exhibit several bright clumps 2 — 3" from the star 
(Fig. 1). The [O III] and Ha images show two well- 
defined bright clumps north and northwest of the star 
(marked A and B, respectively, in Fig. lb), and a more 
extended clump northeast of the star (clump C). This 
extended clump is especially prominent in the [S II] 
images, where clump A cannot be seen at all. 

We discovered the clumpy structure of the S3 core 
in our earlier study in Paper I, based on uncalibrated 
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Figure 1: (a) [O III] image of the nebula surrounding the WO star, adopted from Paper I. The frame 
shows the MPFS field of view for the observations described in this paper, (b) MPFS observations. 
Isophotes show the brightness distribution in the C IV 580lA line emitted by the WO star. Also 
indicated are the positions of emission clumps, (c), (d), (e), and (f): Images and isophotes of the nebula 
in the He II 4686A, [O III] 5007A, Ha, and [S II] 6717+673lA emission lines, respectively. The cross 
indicates the position of the WO star. Surface brightness is in units of 10~^^ erg cm~^ arcsec"^. 



[O III] and Ha images of the region. However, hav- 
ing no emission spectra for the clumps, we mistook 
the WO star for a clump (clump No. 2 in Paper I). 
Clumps A and B from our new observations corre- 
spond to clumps 3 and 1, respectively, from Paper I; 



clump C, which can be confidently identified only in 
the [S II] lines, was not discovered earlier in the [O HI] 
and Ha images. 

We obtained the integrated spectra of clumps A, 
B, and C; the positions of the image elements whose 
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AX (arcsec) AX (arcsec) 

Figure 2: Ionized gas velocity field in the (a) H/? and the (b) [O III] 5007A lines based on the MPFS 
observations. The cross indicates the position of the WO star. The velocity scale is in kms~^. 



spectra were co-added in each clump are shown in 
Fig. lb. Table 3 gives the relative line intensities for 
the three clumps, corrected for selective extinction. 
We assumed c(H/3) = 0.1 and E(B-V)=0.069, in ac- 
cordance with KB95. For comparison, Table 3 also 
presents the data of KB95 for the integrated spec- 
trum of the central region of S3. 

2.4. Ionized Gas Velocities in the 
Core of S3 

Wc used the rcsiilting spectra to construct the radial- 
velocity field for the central part of the nebula S3. 
We determined the velocity fields from measurements 
made in the main spectral lines from Gaussian fits 
to the corresponding line profiles. Analysis of the 
night-sky spectrum showed that, after linearization 
of the spectra, there remained systematic velocity 
variations in the night-sky lines across the field of 
±0.1 — 0.15 channels, i.e., about 15kms~^. We cor- 
rected our emission-line velocity fields for this sys- 
tematic offset; Fig. 2 shows the results. 

We used the resulting velocity fields to estimate 
the mean velocities in various lines for clumps A, B, 
and C and for the nebula region projected onto the 
WO star. In this last case, we used the same set of 
lenses used when co-adding the spectrum of the star 
(see Section 2.2.). The results are summarized in Ta- 
ble 4, where the quoted errors correspond to the ve- 
locity dispersion within individual spatial elements. 

The resulting mean clump velocities are consis- 



tent with the mean velocity of the nebula S3 as a 
whole derived from our Fabry-Perot interferometric 
measurements {Vhel = — 216 ± lkms~^. Paper I, 
taking into account possible systematic absolute cali- 
bration errors of 10 — 20 kms~^ for the interferometric 
observations), and are in full agreement with the mea- 
surements of Tomita et al. (1993), based on a single 
spectrogram of S3: Vhel = —235 to — 230kms~^. 
We show below (see Section 3) that the velocities 
of the densest H I clouds in the vicinity of S3 are 
w — 230kms~^ or differ from this value by less than 
lOkms^i. 

Two cone-shaped regions where the velocities are 
systematically different from the mean are immedi- 
ately apparent in Fig. 2. Our [O III] line measure- 
ments enable identification of the region where the ve- 
locities differ from the mean value by 40 — 50 kms~^. 
Both deviations from the mean velocity agrc;e with 
the velocity variations in the central part of the bright 
nebula S3 presented in Paper I, based on interfero- 
metric Ha measurements. Note that the cone-shaped 
regions originate at the position of the star, a pattern 
that, generally speaking, might indicate an asymmet- 
ric stellar wind. However, firm conclusions require fur- 
ther observational evidence, since such velocity vari- 
ations could also be a consequence of the inhomoge- 
neous structure of the interstellar medium in the near 
environment of the star. 



3. Radio observations in the 21-cm line 
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Table 2: Fluxes and equivalent widths of emission lines in the spectrum of the WO star. 



A 


line 


F {erg cm s ^) 


EW (A) 


(A) 




This paper 


KB95 


This paper 


KB95 


4659 


C IV+He II 


2.56 X 10"" 


1.48 X 10"" 


360 ± 25 


330 ± 20 


5290 


VI 


9.50 X 10^1^ 


1.26 X 10-1'^ 


13 ±2 


55 


5411 


He II+C IV 


1.28 X 10-1^ 


5.5 X IQ-^^ 


17 ±2 


30 


5470 


C IV 


3.46 X 10"^^ 


4.2 X 10"^^ 


4± 1 


20 


5590 


V 


4.04 X 10-^'^ 


2.46 X 10-15 


62± 11 


120 


5801 


C IV 


2.10 X 10"" 


1.41 X 10-1"* 


420 ± 97 


690 ± 40 



Table 3: Mean velocities for clumps A, B, and C and toward the WO star. 



A 


Line 




l/l{Hp) = 


100 




(A) 




A 


B 


C 


KB95 


4686 


He II 


31.1 ±0.7 


28.3 ±0.8 


12.2 ±0.8 


23.1 


4711 


[Ar IV] 


7.1 ±0.7 


7.9 ±0.8 


3.1 ±0.8 


4.5 


4740 


[Ar IV] 


6.3 ±0.7 


6.5 ±0.7 


5.8 ±0.8 


5.62 


4861 




100.0 


100.0 


100.0 


100.0 


4959 


[0 III] 


174 ±5 


168 ±8 


164.0 ±2 


188 


5007 


[0 III] 


519 ± 13 


511 ±20 


489 ±6 




5876 


He I 


4.1 ±0.3 


4.6 ±0.7 


5.4 ±0.8 


7.1 


6542 


[N II] 


1.4 ±0.2 


1.8 ±0.3 


2.3 ±0.3 




6563 


Ha 


273 ± 4.0 


283 ± 10 


273 ±5 


276 


6583 


[N II] 


2.4 ±0.2 


2.8 ±0.4 


4.6 ±0.3 


6.4 


6717 


[S II] 


8.9 ±0.3 


10.2 ±0.2 


15.2 ±0.9 


8.73 


6731 


[S II] 


6.9 ±0.3 


7.6 ±0.2 


10.4 ±0.8 


5.78 


/(i?/?) erg cm ^ s ^ 


1.5 X 10"'* 


1.2 X 10"'* 


7.7 X 10"^'' 





3.1. Observations and Data Reduction 

We constructed a high angular resolution VLA H I 
map covering an extended region of the galaxy IC 
1613; inchiding the WO star and the associated bipo- 
lar shell (this was part of a large project to investi- 
gate the neutral-gas structure and kinematics in IC 
1613 by Wilcots ct al., in press). The maps shown 
below were constructed from combined B, C, and D- 
array VLA observations with a velocity resolution of 
2.57 km s""'^. The data were Hanning smoothed, sub- 
jected to the usual calibration procedure, and made 
into maps using the AIPS package. The results are 
represented in the form of a data cube with angular 
resolution 7. "4 x 7."0, which corresponds to a linear 
resolution of -^23 pc. 

Figure 3 shows the H I distribution for the east- 
ern sector of the galaxy constructed using these data. 
A well-defined giant cavity (supercavity) surround- 
ing a large fraction of the stellar population of IC 
1613 is immediately apparent. This giant, low surface 
brightness region of 21-cm emission is also visible on 



the map of Lake & Skillman (1989). Earlier Paper I, 
we suggested that this region was deficient in neutral 
gas because it had been swept out, implying that the 
supercavity should be surrounded by a dense, wind- 
blown shell. Figure 3 shows that such a shell or ring 
of neutral gas does, indeed, exist, and that the WO 
star resides in a dense gas layer near its inner edge. 

The angular size of the supercavity is about 5', 
corresponding to a linear size of 1 kpc. The char- 
acteristic thickness of the shell is Ai? ~ 1 — 2', or 
200-350 pc. 

A sketch of the bipolar ionized shell surrounding 

the WO star and its relative position with respect to 
the H I supercavity are shown on an enlarged scale in 
Fig. 4, which presents the southeastern sector of the 
H I ring. This figure shows that the core of the ionized 
bipolar shell (previously known as S3) and its bright 
southeastern component are located in a dense region 
of the ring, whereas the extended, faint northwestern 
component is in the low-density medium inside the 
supercavity. 



01 02^ 40^ 30^ 20^ 10^ 

RA (B1950) 

Figure 3: Eastern sector of the galaxy IC 1613: brightness distribution in the 21-cm line from the 
VLA data (grey scale) and Ha isophotes (based on the observations made with the 4-m KPNO 
telescope). An enlarged-scale version of the framed region surrounding the nebula S3 and WO star 
is shown in Fig. 4. 

Table 4: Mean velocities of clumps A, B, C, and toward the WO star. 



Lines 




Vhel, 


kms ^ 






A 


B 


C 


WO 


Hell 4686 


-213 ± 14 


-227 ±22 


-186 ±30 


-214 ±21 




-250 ± 10 


-242 ± 13 


-252 ±7. 


-231 ± 18 


[OIII]4959 


-247 ± 11 


-241 ± 13 


-257 ±6. 


-230 ± 18 


[OIII]5007 


-241 ±8. 


-238 ± 17 


-253 ± 10 


-225 ± 18 


Ha 


-246 ± 6. 


-248 ± 5. 


-242 ± 4. 


-251 ±6. 


[SII] 6717 


-247 ±23 


-225 ± 13 


-224 ± 11 


-253 ±26 


< Vhel > 


-241 ± 6 


-237 ±4 


-236 ± 11 


-234 ±6 
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01 02 40 35 30 25 20 

RA (B1950) 




01" 47' 



01^ 02"^ 32^ 30^ 28^ 26^ 24^ 22^ 

RA (B1950) 

Figure 4: top - Southeastern sector of the dense H I ring. The contours correspond to H I 
column densities of N(H I) = (2,6,18) x lO^o cm . The cross indicates the position of the 
WO star, bottom - Sketch of the bipolar ionized shell associated with the WO star (isophotes 
adopted from Paper I); the Ha isophotes correspond to brightness levels of 0.028, 0.2, 0.4, 1.0, 
and 2.0 x 10~^^ erg cm~^ s~^ arcsec"^. The H I distribution is superimposed as a grey-scale image 
(including the framed region in Fig. 3). 
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Figure 4 clearly shows a local decrease in the 21- 
cm brightness, filled by the ionized gas of the core 
and southeastern component. Wc can also sec two 
prominent, extended H I features just on the outside 
of the two sides of the northwestern component of the 
ionized shell, which reproduce the shape of the shell. 
This structure probably provides evidence for a local 
neutral shell surrounding the bipolar shell. 

The parameters of the most prominent H I clouds 
forming the local neutral shell surrounding the ion- 
ized bipolar feature are summarized in Table 5, which 
presents the central coordinates, total mass, and max- 
imum column density of each cloud. 

3.2. Neutral Gas Velocities 

In order to reveal the possible effect of the expansion 
of the local neutral shell, we constructed the radial- 
velocity distribution within a right-ascension band 
passing through the WO star (Fig. 5). The band is 
200" wide and therefore includes all clouds making 
up the local shell. This distribution clearly demon- 
strates velocity variations with distance from the star 
that are typical of an expanding shell (the so-called 
'velocity ellipse') at right ascensions from 0l''02"'28'' 
to 0l''02'"31^ This region corresponds to the part of 
the local H I shell that is east of the WO star; the 
characteristic radius of this 'half-shell' is 100-140 pc. 
The velocity difference between the two sides of the 
neutral shell is ^ 20kms~"'^, suggesting an expansion 
velocity of ~ lOkms"^. The velocities of the densest 
'unaccelerated' H I clouds in a large region surround- 
ing the WO star are \Hei = -230 ± Skms'^ 

4. Discussion 

The results of spectroscopic observations of the WO 
star and the central region of the nebula S3 reported 
in this paper are in full agreement with earlier ob- 
servations, within the errors, when the different slit 
locations are taken into account. The spectra of the 
compact clumps identified in this paper do not dif- 
fer significantly from the those of the nebula S3 (see 
Garnett et al. (1991), KB95 and references therein). 
The relative line intensities listed in Table 3 are sug- 
gestive of a typical H II region ionized by a star with 
high effective temperature. The observed differences 
between the spectra of the three clumps could be due 
to different excitation conditions; distance from the 
exciting star and electron density. No shock signa- 
tures were found either in the emission spectra (the 
[N II] and [S II] lines are weak) or in the mean veloc- 
ities of the bright clumps. The ionized gas velocities 



-40 -20 20 40 60 




-180 -ZOO -ZZO -240 -Z60 

v(hei) km/s 

Figure 5: HI radial- velocity variations in a right- 
ascension band from the VLA observations (see 
text). The fluxes are in mJy/beam. 
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Table 5: Parameters of HI clouds in the vicinity of the bipolar shell. 



Clouds 


RA(B1950) 


Dec(B1950) 


N(H1), cm-'' 


M, Mq 


a 


1 02 28 


1 47 55 


5 X 10^" 


~ 4 X 10-* 


b 


1 02 31.5 


1 48 20 


8 X 10^0 


~ 3 X 10^ 


c 


1 02 21-24 


1 48 08 


2 X 10^1 


1 X 10-^ 


d 


1 02 26 


1 48 20 


7 X 1020 


2 X lO'' 


e 


1 02 28 


1 48 30-60 


9 X 10^° 


6 X 10^ 



of the clumps in the central core of the nebula agree 
with velocity measurements of neutral gas clouds in 
the vicinity of the bipolar shell. 

4.1. Electron density and He abundance 
in the nebula core. 

The relative intensities of the [S II] doublet can be 
used to estimate the electron density in the three 
bright clumps identified in this paper. Wc determined 
the electron densities in clumps A, B, and C us- 
ing formulas adopted from De Robertis et al.(1987), 
assuming that the electron temperature of the neb- 
ula is Te = 17000 K (as determined in KB95 from 
[O III] 4959/4363 A line intensities). The resulting 
values are listed in Table 6. 

We derived the He abundance in the nebula from 
He I and He II line intensities of individual clumps 
using the diagnostic relations and recombination co- 
efficients for hydrogen and helium from Kholtygin 
(1997),. We determined the relative heliiun abun- 
dances N(He+)/N(H) and N(He++)/N(H) and the 
total helium abundances N(He)/N(H) (in number) of 
the clumps. The results arc given in Table 6, which 
also lists data adopted from KB95. We can see that 
our abundances coincide with those of KB95, and are 
in good agreement with the estimate of Davidsonfe 
Kinman (1982): N(He)/N(H) = 0.073 ± 0.009. (In 
both Davidson & Kinman (1982) and KB95, the He 
abundances were obtained from the spectrum of the 
entire bright core of the nebula without subdividing 
it into individual clumps). 

4.2. Velocity of the stellar wind of the WO 
star. 

One of the major parameters determining the struc- 
ture and kinematics of the interstellar medium at 
large distances from WR stars is the stellar-wind ve- 
locity. Kingsburgh ct al.(1995) have shown (see also 
KB95) that the halfwidth of the C IV 5801 A line 
at zero intensity is a good estimator of the ter- 



minal wind velocities of WR stars. The C IV line 
halfwidth at zero intensity measured from our spec- 
tra for the WO star yields a terminal wind velocity of 
Vyj = 2910kms~^, which is virtually coincident with 
the value 2850 kms~^ obtained in KB95, given that 
the accuracy of our mcasiircmcnts is determined by 
the width of the spectral channel (~ 130kms~^). We 
thus confirm the conclusion of KB95 that the termi- 
nal wind velocity of the WO star in IC 1613 is the 
lowest among the six objects of this class currently 
known in the Local-Group galaxies. 

4.3. Large-scale Structure of the Interstel- 
lar Medium. 

In principle, the giant H I ring shown in Fig. 3 could 
be the projection onto the plane of the sky of a shell or 
torus. A torus could be the result of the wind-blown 
siipcrcavity breaking out of the gaseous disk of the 
galaxy. In both cases, the neutral gas ring surround- 
ing the supercavity displays a characteristic filamen- 
tary structure. The filamentary system making up the 
H I ring bifurcates at its southeastern end. A sys- 
tem of arc-like filaments with characteristic lengths 
of 1 — 1.5', or 200-300 pc, is also apparent. This fila- 
mentary pattern undoubtedly testifies to the action of 
dynamically active processes that are responsible for 
the formation of the supercavity and the surrounding 
dense ring or shell. This large-scale structure appar- 
ently formed as the result of a burst of star forma- 
tion in the region. However, determining the nature 
of this burst is a separate problem, which lies be- 
yond the scope of this paper. Of interest to us here 
is the very existence of such a dense shell or torus 
surrounding the supercavity, since we consider this 
large-scale mass distribution to be responsible for the 
unique bipolar structure of the extended ionized shell 
associated with the WO star. 

Lozinskaya (1997) suggested that the bipolar 
structure of the ionized shell surrounding the WO star 
might be the result of a strong stellar wind break- 
ing out of the dense gas layer. Earlier Paper I, we 
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Table 6: Electron density and He abundances of 
clumps A, B, and C. 





A 


B 


C 


KB95 


Tie cm 


170 


80 


< 30 


100 


N(Hc+)/N(H) 


0.051 


0.057 


0.067 


0.057 


N(He++)/N(H) 


0.028 


0.025 


0.011 


0.025 


N(He)/N(H) 


0.079 


0.082 


0.078 


0.082 



associated this dense layer with a hypothetical wind- 
blown shell, which wc suggested should surround the 
giant H I-deficient region in IC 1613 seen on the map 
of Lake & Skillman(1989). However, only in this pa- 
per have we finally been able to provide some evi- 
dence supporting this hypothesis. The following fac- 
tors form our principal arguments. 

1. One piece of evidence is the relative orientation of the 
bipolar ionized shell with respect to the H I supercav- 
ity. It is clear from Fig. 4 that the bright core of the 
nebula is extended along the boundary of the dense 
H I ring; the smaller southeastern component of the 
ionized bipolar shell lies inside a dense neutral gas 
layer, and the extended northwestern component lies 
in the low-density medium inside the supercavity. 

2. Afanas'cv ct al. (2000) showed that the southeastern 
component is brighter than the northwestern compo- 
nent, and also has a lower expansion velocity. 

3. The local cavity and extended H I features whose 

discovery is reported in this paper testify to the pres- 
ence of a neutral shell surrounding the ionized bipolar 
shell, providing direct evidence for a physical associ- 
ation of the WO star and the ionized bipolar shell 
with the supercavity and dense H I ring (the only 
alternative being a chance projection). 

The gas density in the clouds making up the local 

H I shell surrounding the bright central and south- 
eastern components of the ionized nebula can be de- 
rived from our maximum column density and total 
mass estimates, assuming that the clouds have the 
same sizes in the line-of-sight and transverse direc- 
tions. The clouds that make up the local neutral H I 
shell have irregular structures. Clouds a, b, and c in 
Table 5, which make up the expanding part of the 
H I half-shell, have lengths in the plane of the sky 
ranging from 40 to 120 pc. For a column density of 
(5 — 9) X 10^" cm~^, this yields a characteristic density 
of njj J ~ 1.5 — 7.5 cm~^. 

For comparison, wc present estimates of the av- 
erage density of ionized gas in this nebula: Davidson 
& Kinman (1982) found Ug = 8.5 cm~^ based on the 



H/3 emission of the central 6" of the nebula S3; Keni- 
cutt (1984) reports rie — 1 cm~^ based on a study of 
the integrated Ha emission of the bright nebula; and 
Goss & Lozinskaya (1995) found Ug = 3.5 cm""^ from 
their analysis of its thermal radio emission. 

If the H I clouds surrounding the bright central 
and southeastern components of the bipolar ionized 
shell are, indeed, parts of a shell of neutral gas blown 
out by the stellar wind, the radius and expansion ve- 
locity of the expanding section east of the WO star 
should be 100 — 140 pc and about 10 kms~^, re- 
spectively (see Section 3.2). In the standard model 
(Castro et al.,1975; Weaver et al.,1975), we find the 
kinematic age of the wind-blown neutral shell to be 
t = (6 — 8) X 10^ yr. (Note that this estimate depends 
neither on the stellar wind power nor on the density 
of the surrounding gas). This age substantially ex- 
ceeds the duration of the final WO stage, and corre- 
sponds to the main-sequence lifetime of a massive star 
(M > 30 - 4OM0). We thus conclude that the local 
nciitral shell surrounding the bright southeastern and 
central components of the ionized nebula could have 
been produced by the stellar wind of the WR pro- 
genitor when it was still a main-sequence star. Note 
that, earlier Paper I, we concluded that the ionized 
bipolar shell was blown out by the successive action 
of the stellar wind from the central star, first at the 
WR and then at the WO stage. 

Our estimates indicate that the terminal velocity 
of the stellar wind of the WO star reaches 3000 kms""'^ 
(see Section 4.2); a mass-loss rate of 2.9 x lO~^M0/yr 
was derived in KB95 based on the relation between 
the He II line luminosities and mass-loss rates for two 
WO stars. Accordingly, the power of the stellar wind 
of the WO star in IC 1613 is Ly^ ~ 10^^ erg/s. In all 
likelihood, the wind power was a factor of two to three 
lower during the earlier WR stage (WR stars have 
typically terminal velocities of 2000 km s~^), but the 
duration of the WR stage (5 x 10''' yr) is at least an 
order of magnitude longer than that of the WO stage. 
The power of the stellar wind of the main-sequence 
progenitor could have been two orders of magnitude 
lower than at the final Wolf-Rayet stage, however, 
the main-sequence lifetime was at least an order of 
magnitude longer than the Wolf-Rayet stage. 

The formation of the wind-blown shell is deter- 
mined by the standard relation: 

R = 2i{Li^/noY^Hl'^ 

, where L36 is the power of the stellar wind in units 
of 10'^^ erg/s, rio the density of the surrounding gas, 
R the radius in pc, and the time in units of 10^ yr 
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(Castro et al.,1975; Weaver et al.,1975). We can eas- 
ily see that, although the mechanical energy inputs 
from the massive star during the WR and WO stages 
are comparable, the stellar wind produced during the 
main-sequence stage is capable of creating a larger 
shell than the powerful but short-lived stellar wind 
at the final stage. Therefore, this latter wind serves 
as an additional source of mechanical energy in an 
already formed, slowly expanding shell. 

Our results for the two extended arc-like H I 
clouds surrounding the northwestern component of 
the ionized shell arc much less conclusive. On the 
one hand, these formations could be the result of the 
stellar wind of the WO star and its progenitor break- 
ing into the tenuous gas inside the supcrcavity. Our 
measurements, indeed, indicate that the H I surface 
density inside the supercavity is at least an order of 
magnitude lower than in the clouds listed in Table 5. 
The thickness of the emitting gas layer can be taken 
to be twice that of the supersheU in the plane of the 
sky; i.e., about 3 — 4' (600-700 pc), according to Fig. 
3. Accordingly, we find the mean gas density in the 
supercavity to be ~ 0.04 — 0.05 cm~^. If the observed 
ring-shaped H I distribution is the projection of a 
thick torus rather than a shell, the size of the emis- 
sion region should be even larger, so that the above 
estimate should be considered an upper limit to the 
gas density inside the supercavity. The same standard 
relation yields a radius of i? 180 pc for a shell blown 
out in the tenuous medium of the estimated density 
over the time t = 7 x 10^ yr found above. This agrees 
with the observed size of the northwestern part of the 
local shell outlind by two extended arc-like H I clouds 
in Fig. 4. 

On the other hand, the formation of such an arc- 
like structure outlined by these extended H I clouds 
could be the result of the collective stellar winds (and 
possible supernova explosions) from association No. 9 
of the list of Hodge (1978), which is adjacent to and 
inside the supercavity. (This association is the largest 
in size and the richest in the galaxy) . In this case, the 
stellar wind and ionizing radiation of the WO star act 
inside this already formed arclikc feature. 

We plan a detailed analysis of the kinematics of 
the ionized and neutral gas components in this re- 
gion, and a comparison with the stellar population. 
This study should lead to unambiguous conclusions 
about the nature of the arclike formation and the 
possible genetic relation between the WO star and 
the association. 



5. Conclusions 

We have observed the central part of the nebula S3 
surrounding the WO star (field spectroscopy with 
the 6-m telescope) and an extended surrounding re- 
gion more than one kpc in size (21-cm VLA ob- 
servations). We have mapped the brightness of the 
central nebula in the main spectral lines and con- 
structed the radial-velocity distribution based on H/3 
and [O III] 5007A line measurements. 

We identified three compact clumps, and deter- 
mined their integrated spectra, mean velocities, and 
electron densities. Our spectroscopic data for both 
the WO star and the central region of S3 are in 
agreement with previous observations (see Garnett et 
al.(1991), KB95 and references therein). The relative 
line intensities in the spectrum of the ncbida arc in- 
dicative of pure recombination emission from a highly 
excitated H II region. 

We analyzed the structure and kinematics of the 
neutral gas in an extended region aroimd IC 1613 
that includes the WO star and the associated ionized 
bipolar shell. A giant, dense shell or toroidal forma- 
tion of neutral gas about 1.5 kpc in size surrounds the 
H I supercavity. This giant structure was very likely 
produced by a burst of star formation in this region. 

We also discovered for the first time a local de- 
ficiency in the brightness of the 21-cm line and two 
extended H I features, providing evidence that the ex- 
tended ionized bipolar shell associated with the WO 
star (Lozinskaya, 1997 and Paper I) is, in turn, sur- 
rounded by a local shell of neutral gas. The mass and 
expansion velocity of this shell indicate that it was 
blown out by the stellar wind of the WO progenitor 
when it was still a main-sequence star. 

Our new observations support a scenario in which 
the extended bipolar structure was produced by the 
stellar wind of the WO star residing in the dense 'wall' 
of the supercavity. The arclike structure of the ex- 
tended H I clouds just outside the faint northwest- 
ern component of the ionized shell might have been 
shaped by the collective stellar wind (and supernova 
explosions?) of Hodge's (1978) association no. 9. In 
this case, the wind and ionizing radiation of the WO 
star act inside this previously formed arclike struc- 
ture. 
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